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ABSTRACT 

We present Spitzer Space Telescope IRAC photometry of the Local Group 
dwarf irregular galaxy IC 1613. We compare our 3.6, 4.5, 5.8, and 8.0 /im pho- 
tometry with broadband optical photometry and find that the optical data do 
not detect 43% and misidentify an additional 11% of the total AGB population, 
likely because of extinction caused by circumstellar material. Further, we find 
that a narrowband optical carbon star study of IC 1613 detects 50% of the total 
AGB population and only considers 18% of this population in calculating the 
carbon to M-type AGB ratio. We derive an integrated mass-loss rate from the 
AGB stars of 0.2-1.0x10"^ Mq yr~^ and find that the distribution of bolomet- 
ric luminosities and mass-loss rates are consistent with those for other nearby 
metal-poor galaxies. Both the optical completeness fractions and mass-loss rates 
in IC 1613 are very similar to those in the Local Group dwarf irregular, WLM, 
which is expected given their similar characteristics and evolutionary histories. 

Subject headings: stars: AGB - stars: carbon - stars: mass loss - galaxies : dwarf 
- galaxies : irregular - galaxies : Local Group - galaxies : individual(IC 1613) 



1. Introduction 



IC 1613 is an isolated member of the Loc al Group and serves a s the prototype fo r the 
DDO Ir V type dwarf irregular (dl) galaxy (Ivan den Berghl 2000al) . Discovered by IWoli 
( 1l906h . IC 1613's relative prox imity (730 kpc; iDolphin et al.ll200lh . low inclination angle 
(i = 38° lLake fc Skillmanlll989l ). favorable galactic latitude (b = —60°), and low foreground 
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reddening (E(B— V)= 0.025 mag; ISchlegel et al.l Il998l ) make it an ideal target for stellar 
population and interstellar medium (ISM) studies. Thi s is evidenced by the impressive 
number of previous studies of this galaxy, summarized by Ivan den Berghl (j2000al ). 



Lake fc SkillmanI (119891 ) studied the neutral hydrogen in IC 1613. They found the galaxy 



to be gas ric h with the H I dis tribu tion being lumpy in th e inner galaxy and smoother at 
larger radii. ICole et al.l (119991 ) and ISkillman et al.l (120031 ) performed deep photometry of 



single HST/WFPC2 fields in the inner and outer galaxy, respectively. Both of these studies 
used high-precision photometry to derive star formation histories. They found continuous 
star formation over the past few hundred Myr, as well as populations of intermediate age 
and ancient stars. 

IC 1613 is, in most respects, a typical irregular galaxy. Its metal abundance, current 
star formation rate, and dust content are all normal for a galaxy of its luminosity. The 
morphology consists of a prominent star forming region northeast of the galaxy center and a 
bar running linearly through the galaxy center from southeast to northwest. P erhaps the only 
exceptional property of I C 1613 is its nearly complete absence of star clusters ( van den Bergh 



19791 



2000bl ) 



Wyder et al.ll2000t ). which may be evidence of its isolated evolution ( 



van den Bergh 



In this paper we present Part II of our Spitzer IRAC thermal infrared study of the 
stellar populations in Local Group dwarf irregulars with a census of the asymptotic giant 
branch stars in IC 1613. In §5] we describe the observations, data reduction, and photometry. 
In §3] we present the Spitzer IRAC images and describe the near- and mid-IR morphology 
of IC 1613. We characterize the infrared color-magnitude diagram (CMD) in [|H In §5] 
we detail a complete census of the asymptotic giant branch (AGB) stars, a comparison with 
broadband optical photometry ( §5.11) . a comparison with an optical carbon star study ( §5.2p . 
the determination of mass- loss rates (MLRs; §5.3p . and the total luminosity contribution from 
the AGB population ( §5.41) . In §6] we summarize our conclusions. 



2. Observations and Data Reduction 



2.1. Infrared and Optical Data and Photometry 



Images of IC 1613 were obtained with t he Infrared Array C amera (IRAC; iFazio et al. 



20041 ) on board the Spitzer Space Telescope ( IWerner et al.l 120041 ) in two separate AORs on 
2003 December 21 UT and 2005 December 23 UT (AORKEYs: r5051648 and rl5892736), 
respectively. Both AORs are part of a larger guaranteed time observing program (PID: 128, 
PI: R.D. Gehrz). Five individual 200 s frames at 4.5 and 8.0 fim and 32 individual 30 s frames 
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at 3.6 and 5.8 yum were used resulting in total integration times of 1000 s per pointing at 4.5 
and 8.0 //m and 960 s per pointing at 3.6 and 5.8 fim. These data yielded 5-a point source 
sensitivities of 4.5, 3.5, 10.7, and 16.5 /iJy at 3.6, 4.5, 5.8, and 8.0 fim, respectively. Because 
the IRAC 3.6 and 5.8 /im field of view is offset from that at 4.5 and 8.0 fim, a 3x2 mosaic 
in array coordinates approximately 15'xl0' for each channel was created to map IC 1613. 
This resulted in a coverage area with data in all four IRAC channels of approximately 9'x9' 
in size, centered at a(2000)=02:06:30, (5(2000)=01:04:50. Figure □ shows this field of view 
overlaid on the Digitized Sky Survey image of IC 1613. 

IRAC data (pipeline version 13.2.0) were reduced with the MOPEX0 reduction pack- 
age, version 2006 March 1. The backgrounds of individual frames were matched using the 
overlap routine. The mosaic program was used for outlier detection, image interpolation 
and co-addition, and to create the final mosaics. Very bright objects imaged in the program 
immediately preceding our observations resulted in persistent images in the 8.0 fim data 
taken in 2003. These artifacts were effectively removed by subtracting a median-combined 
image of all of the dithered frames from each individual frame. 



Point source photometry was done using the DAOPHOT II photometry package (jStetson 



19871 ). A sharpness clip was applied to the photometry in each of the IRAC bands. This clip 
removes any objects from the photometry list whose radial profiles are significantly steeper or 
broader than the profile of the derived point spread function (PSF). We adopt the standard 
DAOPHOT errors as the total photometric errors for our sources. T he IRAC instrument 



also has an estimated absolute calibration accuracy of 3% (iReach et al.ll2005l ). although we 
do not include this uncertainty in our photometric errors. 

Op tical V and / photom etry of IC 1613 were obtained from the OGLE-II microlensing 



survey (lUdalski et al. 



20011 ). These data are used to aid in identifying the stellar type of 



the stars detected in our IRAC data, which is challenging with IRAC data alone (see §1]). 
The optical data were acquired in the context of detecting objects that are variable (due to 
either microlensing events or intrinsic stellar variability) by taking a large number of frames, 
each with relatively short integration times. The photometry presented here is the mean of 
the individual phot ometric observations , including those objects found to be variable. The 
field of view of the lUdalski et al.l (120011 ) optical data covers the entire IRAC 3.6, 4.5, 5.8, 
and 8.0 fim coverage area, except for a small triangular region ~0.9 square arcmin in size in 
the south-eastern edge of the IRAC coverage area (see Figure [1]). 



^MOPEX is available from the Spitzer Science Center at 



http : //ssc . spitzer . caltech.edu /postbcd/ 
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2.2. Foreground Star Contamination 



To estimate the number of foreground stars i n both the optical and IR data we used the 
Milky Way stellar population synthesis model of iRobin et al.l ( 120031 ) . Sources were modeled 
for t he Galactic coordin ates of IC 1613 in a one square degree field to reduce statistical error. 
The iRobin et al.l (120031 ) model provides magnitudes in the Johnson-Cousins L-band, whose 
central wavelength is close to that of the IRAC 3.6 /xm band. In the region of IC 1613 where 
we have coverage in the optical and all four IRAC channels, we expect 15 foreground stars 
in our IR color- magnitude diagram (CMD). These 15 foreground stars are expected to have 



-7.5<M3.6^„ < 



-16 and [3.6] — [4.5] colors near zero. 



3. Infrared Galaxy Morphology 



IRAC 3.6, 4.5, 5.8, and 8.0 /xm images of IC 1613 are shown in Figure [2l In general, 
the stars in these images are uniformly distributed througho ut the images, consistent wit h 
the regular elliptical distribution of red giant s discuss e d bv iBaade fc GaposchkinI (119631 ) . 
and dubbed Baade's Population II 'sheet' by ISandagd (Il97ll ). In addition to the smooth 
distribution of stars, a concentration of bright stars is seen in the large star forming complex 
in the northeast. The only diffuse emission detected in these images is 8.0 /im emission 
from hot (~350 K) dust and/or polycyclic aromatic hydr ocarbons. This emis sion is detected 
in regions roughly coincident with strong Hck emission (j Jackson et al.ll2006l . and references 
therein) . 



4. Comparison of Optical and IR Photometry 



The optical color-magnitude diagram (CMD) from the lUdalski et al.l (120011 ) data is 
shown in Figure [3l segregated into regions consisting of (a) blue objects that can include 
massive main sequence stars and unresolved H II regions, (b) asymptotic giant branch (AGB) 
stars, (c) red supergiants (RSGs), and (d) red giants (RGs) below the tip of the red giant 
branch (TRGB). We refer hereafter to the stars called out in (d) as sub-TRGB stars. These 
regions are used to identify the stellar populations in the infrared CMD, since IRAC data 
alone provide little information on stellar effective temperatures. These indicated regions are 
separated by gaps reflecting the estimated 1-cr photometric error, which were employed to 
minimize misidentification of stellar types solely due to photometric errors. Reddening from 
circumstellar and interstellar material can also cause such misidentification. The optical 
CMD shows features of both recent and very old star formation, consistent with the finding 
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of ISkillman et al.l (120031 ) that IC 1613 has a relatively constant star formation rate over the 
age of the Universe. 

Figure m shows the 3.6 /im absolute magnitude versus [3.6] — [4.5] fim CMD. The left 
panel shows all of the stars detected at both 3.6 and 4.5 /xm, with photometric errors, 
averaged over one magnitude bins, displayed. The left panel also shows the 50% completeness 
limit (black line) derived from artificial star tests. The right panel shows the stars detected 
in V, I, and at 3.6 and 4.5 /um, with RSGs shown as black stars, AGB stars as green triangles, 
red giants as red squares, and blue objects as blue circles. The right panel contains a vector 
showing the dis placement a star would experience due to 10 vis ual magnitudes of extinction 
and reddening (llndebetouw et al.ll2005l : iRieke fc Lebofskylll985l ). though this vector does not 
include the effects of dust emission, which can be significant in the mid-IR. The rig ht panel 
also shows the AGB limit, which was determined by taking the iGroenewegenI (120061 ) models 
for carbon-rich AGB stars with Teff = 2650 and Tg// = 3600 and negligible mass- loss and 
scahng them to an absolute bolometric magnitude of M^oi = —7.1. The AGB limit depicted 
in this figure is a line connecting these two models. AGB stars that are losing significant mass 
can be brighter than this limit due to thermal emission of the mass-loss material increasing 
their 3.6 //m flux (as shown in detail in §5.3p . however, any such objects would also have 
red [3.6] — [4.5] colors. None of these objects are observed, thus, we conclude all of the stars 
brighter than the AGB limit are true red supergiants. 

The right panel also shows the 3.6 fim TRGB. The value of the TRGB we adopt 
(M3 6=— 6.2±0.2) was determined by inspecting the 3.6 fim luminosity function for sources 
detected at both 3.6 and 4.5 /im. The 3.6 /im luminosity function (the number of stars 
in each 0.2 magnitude bin at 3.6 /zm) is shown in Figure [51 We adopt this value of the 
TRGB based on the abrupt drop in detections at that magnitude. Also, optically classified 
sub-TRGB red giants with blue [3.6] — [4.5] colors are observed with M3.6 magnitudes up to 
this value, but above this value only optically classified sub-TRGB red giants with very red 
colors are observed. The IR fluxes of these objects are consistent with mass-losing AGBs 
(see ^ rather than sub-TRGB red giants, supporting our adopted value of M3.6=— 6.2 for 
the 3.6 fim TRGB. This value is 0.2 r aagnitudes fainter tha n the value found for the Large 
Magellanic Cloud (LMC ) (Mf/ =— 6.4: Ivan Loon et al.ll2005l ) and 0.4 magnitudes lower than 
that adopted for WLM (j Jackson et al.l 120071 ). It is unlikely that the difference in values of 
the TRGB between IC 1613 and WLM are actually so disparate, given their similar metal- 
licities. The difference in these values likely reflects the uncertainty in our adopted TRGB 
values. This uncertainty is not a major concern, however, because shifting the value of the 
TRGB by up to 0.5 magnitudes affects the detection statistics only slightly (see §5.11) . 

In general, the M3.6 versus [3.6] — [4.5] CMD of IC 1613 is very similar to that of WLM, 
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with the only major difference being that there are significantly more sub-TRGB red giants 
detected in IC 1613 than in WLM because WLM is more distant. As is the case of WLM, a 
small population of stars redward of the main stellar distribution is detected, though fewer 
are seen in IC 1613 than in WLM. As we discuss in §5.31 the infrared fluxes and colors of 
these objects are consistent with mass-losing AGB stars. 

Figure lis the Mg.o versus [3.6]-[8.0] CMD for IC 1613. This CMD shows a nar- 
row vertical feature with [3.6] — [8.0] ~0 and another broad distribution of red objects with 
—12 < Mg o < —7.5 1 < [3.6] — [8.0] < 4. These red objects are also the most luminous, 
reddest objects in the M3.6 versus [3.6] — [4.5] CMD, with typical colors between 0.5 and 1.0. 

Figure [7] shows the spatial distributions of different stellar types based on their optical 
and IR fluxes. The distributions of both red giants and AGB stars are very smooth with 
no obvious concentrations aside from the radial stellar gradient. A thin gap is observed in 
the optical da t a, wh ich is certainly an instrumental effect, though this is not mentioned in 



Udalski et al.l (I2OOII ). There are two conspicuous features in the blue objects, which outline 



the young stellar distribution; a bar running across the midplane of the galaxy and a large 
clump in the northeast, which highlights the most active star forming complex in IC 1613. 
It is clear from this figure that neither the optical nor the IR data presented here reach the 
outermost regions of the stellar population, as the stellar density is still dropping off at the 
edges of the images. In the analogous plot in IJackson et al.l (120071 ). the effects of crowding 
could clearly be seen in WLM as clumps of bright, young stars coincided spatially with an 
absence of fainter red giants. In IC 1613 we se e no evidence for such a n effect, likely because 
I C 1613 is seen mu ch more face on {i = 38°, iLake fc SkillmanI 119891 ) than WLM {i 



69° 



Jackson et al.l 120041 ). 



In Table [2] we present the detection statistics for the optical and infrared photometry 
of IC 1613. We dete ct 5211, 3162, 8 86, a nd 618 point sources at 3.6, 4.5, 5.8, and 8.0 fim, 
respectively, and the lUdalski et al.l (1200 ll ) V and / data contain 7266 point sources within 
the optical and IR coverage region. We detect 183 point sources in all four IRAC channels 
with no corresponding detection in the optical. Three of these are RSGs (i.e., they are above 
the AGB limit), two of these are very red objects below the TRGB, and the remaining 178 
are either RSGs, AGB stars, or blue objects. 



5. The AGB Stars 



As shown in Figure Hj, our completeness limit is more than 0.5 mag fainter than the 
TRGB for the bluest objects and 1.5 mag fainter than the TRGB for red objects. Con- 
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sequently, our data constitutes a complete census of the super-TRGB AGB stars with the 
exception of a very small number of objects not detected due to crowding. 



As described in IJackson et al.l (120071 ) . though our data should represent a complete 
census of the AGB stars, we cannot distinguish between AGB stars, RSGs, and blue objects 
using IRAC data alone. However, above the TRGB, optically identified RSGs and blue 
objects comprise only 9% (95 out of 1041) of the objects between the TRGB and the AGB 
limit. Therefore, we assume any object between the TRGB and the AGB limit is an AGB 
star unless its optical identification dictates otherwise, with the understanding that there 
may be a small amount of contamination by other stellar types. 



5.1. Optical Completeness 



Altogether w e dete ct 1052 AGB stars in our IRAC data. Of these 57% are detected in 
the lUdalski et al.l (120011 ) V and / data. In addition to the 43% of AGB stars detected in the 
IR but not seen in the optical, 11% of the IR detected AGB stars are located in the optical 
CMD where they would be misidentified as sub-TRGB red giants. 



For comparison, in WLM IJackson et al.l (120071) found 3 9% of the IR detected AGB 
stars were not seen in the optical data of iMassey et al.l (120061 ) . and an additional 4% were 
optically detected, but misidentified as sub-TRGB red giants. An i mpor tant check of our 
AGB statistics for WLM was recently provided by IValcheva et al.l (120071 ) . who performed 
J and Ks near-IR photometry. Altogether they detect 355 AGB stars (carbon + M-type 
AGBs), in comparison with the 331 AGB stars we detect in the same field of view. It is 
not surprising that the number of AGBs we detect in WLM and IC 1613 are in such close 
agreement given the striking similarities in the evolution, and particu larly in the intermed iate 
age star formation histories, of these galaxies (Figures 25 and 35 in lDolphin et al.ll2005l ). 



In Figure El we show the fraction of optical completeness as a function of [3.6] — [4.5] 
color (left panel) and as a function of the modeled AGB wind optical depth (right panel; 
AGB mass-loss modeling is described in detail in §5.3p . This figure shows the optical com- 
pleteness fraction is a smoothly declining function of [3.6] — [4.5] color, as is expected if the 
incompleteness is due to extinction by circumstellar material. Much like in the case of WLM, 
we see the completeness fraction fall to 0% at a [3.6] — [4.5] color of ~0.9. 

It is important to understand to what degree our adopted value of the TRGB affects 
these completeness statistics. By shifting our adopted value 0.5 mag fainter the fraction 
of optically undetected AGB stars changes by only 2%, while shifting the TRGB 0.5 mag 
brighter would increase the optically undetected AGB fraction by 10%. As in the case 
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of WLM, this increase in optical incompleteness with increasing luminosity supports our 
conclusion that the optically undetected AGB stars are not seen because they are enshrouded 
by circumstellar material, since the mass-loss phase in AGB stars occurs somewhat above 
the TRGB and should consequently be more conspicuous at higher luminosities. 



5.2. The Carbon Stars 



Narrowband a nd broadband optica l photometry of IC 1613 using the CN/TiO technique 
was perforrn e d by lAlbert et al.l (120001 ) . In comparison with our IRAC photometry, the 
Albert et al.l (120001 ) data detected 50% of the total AGB population seen in our IRAC data. 
This value is slightly lower than the fraction detected by lUdalski et al.l (120011 ) and is caused 
by the brighter com pleteness limits of the narrowband data, si nce the broadband d ata from 
Albert et al.l (120001 ) is significantly more sensitive than that of lUdalski et al.l (120011 ). 



Albert et al.l (120001 ) distinguish between oxygen- and carbon-rich AGBs using a color- 



color diagram. At relatively red (R-I) colors oxygen- and carbon-rich AGBs have easily 
distinguishable CN-TiO colors, however, AGBs with bluer (R-I) colors are not easily sep- 
arated. Consequently, a color limit is imposed, redward of which t he discrimina t ion o f 
chemical composition can be made. Within our IRAC coverage area lAlbert et al.l (|2000l ) 
identified 74 carbon stars. Assuming their C/M ratio of 0.64, this yields 190 total AGB 
stars (C + M-type). Thus, after applying their imposed color limit, 18% of the total AGB 
population was considered for the C/M ratio. 



Figure M shows our IRAC photometry of the carbon stars identified by lAlbert et al. 



( I2OOOI ). The carbon stars have relatively blue colors, probably due to the low-metallicity, 
and consequently the low optical depth of circumstellar wind material of AGBs in IC 1613, 
though four red carbon stars a re seen that are like ly losing significant mass. There are also 
three carbon stars identified by lAlbert et al.l (I2OOOI ) that we find to be well below the TRGB 
at 3.6 /im. These three stars are also below the TRGB in the optical, and are possibly 
extrinsic carbon stars. Because it is very unlikely that these three objects are in fact true 
carbon-rich AGB stars we do not consider them in any of the car bon star stat i stics. With 
the exception of these three objects, the carbon stars identified by lAlbert et al.l (I2OOOI ) have 
relatively bright 3.6 fim luminosities compared to the rest of the AGB stars. 

The comparison between our IR photometry and the narrowband optical data in IC 1613 
is nea rly identical to the values we found for this same c omparison in WLM ( 



20071 ). In WLM we found the completeness fraction of the lBattinelli fc Demerd (|2004l ) optical 



ackson et al. 



narrowband data was 63% of the total AGB population observed in our IRAC data and the 
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C/M ratio was based on 18% of the IRAC detected AGB stars. Like in the comparison of 
broadband optical completeness, the striking similarity between AGB populations in IC 1613 
and WLM is likely due to their similar star formation histories and metallicity evolution. 



In IJackson et al.l (120071 ) we stated that the iBattinelh k Demerd (j200J) optical CN/TiO 
study of WLM detected only 18% of the total AGB population, which is misleading. Again, 
because of the necessity to impose a color limit to discriminate between carbon- and oxygen- 
rich AGBs, relatively blue AGB stars are ignored in the calculation of the C/M ratio. In 
WLM 18% of the total AGB population was used to calculate the C/M ratio, while an ad- 
ditional 45% of the total population was detected but not used in the C/M ratio because 
the chemical composition could not be determined. This point is critical to consider when 
comparing C/M ratio s obtained using diffe rent methods. As an example, recent near-IR 
photometry of WLM (jValcheva et al.l 120071 ) found a C/M ratio of 0.58 (corrected fo r fore- 
ground contamination). This is in sharp contrast with the iBattinelli fc Demerd (120041 ) value 
of 12.4, though these value s must not be directly compared because, as is shown here, the 
Battinelli fc Demerd (120041) C/M ratio was calculated using 18% of the AGB stars, while the 
( IValcheva et al.l 120071 ) value considered the entire AGB population. 



5.3. Mass-Loss Rates 



The mass-loss rates (MLRs) of AGB stars in IC 16 13 were determined b y comparing their 
infrared fluxes with the radiative transfer models of iGroenewegenI ( l2006l ). The [3.6] — [4.5] 
color of AGB stars corresponds to a unique wind optical depth (r), given assumptions about 
the wind composition and stellar effective temper ature {Tpff), as we discuss below. The 
[3.6] — [4.5] colors were linearly interpolated ont o the IGroenewegenI (120061 ) models to estimate 

1 3 

r, then scaled according to the Ivan LoonI (120001 ) MLR prescription, where: M oc r'0~2L4, tp is 
the dust-to-gas ratio, and L is the stel lar luminosity. F or IC 1613 we adopt the nebular oxygen 
abundance of 12+log(0/H) = 7.62 (ILee et al.l 120031 ) . which corresponds to a metallicity of 
0.9% Zf7^. Given tha t the dust-to-gas ratio scales as -0 = t/'qIO"!^^/'^] and ipQ = 0.005 
( van Loon et al. 20051 ). we adopt ipic leia = 4.8x10"^. 



It is more appropriate to apply the metallicity of IC 1613 as it was a few Gyr ago when 
the current AGB population was forming, rather t han the current meta llicity. However, 
determining this value empirically is difficult, and ISkillman et al.l (120031 ) find very little 
metallicity evolution in IC 1613 over the past few Gyr, so applying the current metallicity 
to the AGB stars should be reasonably accurate. 



There are two major uncertainties in deriving MLRs for AGB stars using only our 3.6 
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and 4.5 micron data. The first is that there is a degeneracy between r and Tg//, such that 
cool AGB stars losing no mass can have the same [3.6] — [4.5] color as warmer AGBs that are 
losing mass. Because of this degeneracy we are forced to derive MLRs by adopting a single 
Te// for all of the AGB stars, which is certainly not physical. As a consequence, the MLRs we 
find for indiv idual blue AGB stars (i.e., AGB stars with [3.6] — [4.5] < 0.5) should not be over- 



interpreted. iMarengo et al.l (120061 ) showed that atmospheric molecular absorption features 
can also effect the IRAC colors of AGB stars, introducing additional uncertainty to these 
m easurements, a l thoug h these effects have not yet been quantified. However, as emphasized 



m l Jackson et al.l (120071 ). because the majority of the integrated mass-loss for IC 1613 comes 
from sources much redder than [3.6] — [4.5] > 0.5, where the MLR we calculate is much less 
sensitive to T^ff, the integrated MLR for the galaxy as well as those found for red AGB 
stars are robust. 

The second major uncertainty in our MLR calculation is that we have no a priori in- 
formation about the chemical composition of the AGB winds, and therefore, must assume 
a single chemical composition for the entire AGB population. This assumption is also un- 
physical, and the disparity between the integrated MLRs we measure for all of the modeled 
wind compositions and effective temperatures differ by a factor of four. 

Figure [10] shows the MLR versus bolometric luminosity for each AGB and RSG in 
IC 1613 assuming a wind composition of 85% amorphos carbon (AMC) + 15% SiC and 
effective temperatures of Tg// = 2650 K {filled circles) and Tg// = 3600 K {open circles). 
Again, we consider any object detected at 3.6 and 4.5 /xm that is above the TRGB at 3.6 
fim to be an AGB star, except for those optically classified as blue objects, since unresolved 
H II regions lie in the same part of the IR CM P as mass-losing AGBs. In comparison with 
the tracks of increasing wind optical depth (see lGroenewegenll2006l : I Jackson et al.ll2007l ). this 
plot shows that the most heavily mass-losing stars are among the most luminous AGBs. The 
maximum ML Rs observed in IC 161 3 are in very goo d agreement with the maximum MLRs 



seen in WLM ( iJackson et al.l 120071 ) and the LMC (Ivan Loon et al 



Jura 



1999). and somewhat 



(119841 ). We also see a 



higher than the classic single-scattering MLR limit predicted by 
number of RSGs above the AGB limit, none of which show infrared excess. 

The total MLRs we derive for IC 1613 are 2.4x10"^ M© yr'^ for T^ff = 2650 K and 
85% AMC + 15% SiC, 4.3x10"^ for T^jf = 3600 K and 85% AMC + 15% SiC, 1.0 xlO-^ 
Mq yr-i for Tg// = 2500 K and 60% Silicate + 40% Aluminum Oxide, and 8.8 xlO"^ for 
Teff = 2500 and 100% Silicate. These values are approximately a factor of two lower than 
those found for WLM. Considering that a significant fraction of the mass-loss is coming from 
a small number of objects with short mass-loss timescales, the integrated MLR of a relatively 
low-mass dwarf like IC 1613 could be highly time-variable and the consistency between the 



MLRs of WLM and IC 1613 is notable. 



5.4. Luminosity Contribution Prom AGB Stars 



For point sources detected at 3.6 and 4.5 /im we find a total 3.6 /im flux of 72.4 mJy. 
85% of this total flux is from AGB stars and RSGs, which is very similar to the case of WLM, 
where we found 79% of the flux in point sources was from sources brighter than the TRGB. 
The total flux in poi nt sources we me asure for IC 1613 constitutes 79% of the integrated 4.5 
yum flux detected by iLee et al.l (120061 ) . 



Following Ivan Loon et al.l (120051 ) and adopting an age of the current super-TRGB stellar 
population of 2 Gyr we find a total stellar mass of 1.7x10^ M(d. This value is a factor of 



2.6 hig her than the 6.6x10^ Mq estimated by lLee et al.l (120061 ). Because the Ivan Loon et al. 



(120051 ) method of measuring stellar mass assumes a single-aged stellar population and is 
dependent on the age of that population, our stellar mass may be over-estimated since a 
number of stars above the TRGB in IC 1613 are likely RSGs younger than 2 Gyr old. 



6. Summary of Results and Conclusions 



We present Spitzer Space Te/escoj^e thermal near- IR photometry of IC 1613. These data 
are compared with the optical Fand / photometry of lUdalski et al.l (120011 ) and we find that 
43% of the IR detected AGB stars are not detected in the optical and an additional 11% of the 
AGBs are misidentified as sub-TRGB red giants. We show that the optical incompleteness 
fraction is very well correlated with the [3.6] — [4.5] color, which indicates that the optically 
undetected AGB stars have been reddened by circumstellar material beyond the optical 
completeness limits. 

Our IR p hotom etry is also compared with the narrowband optical carbon star study of 
Albert et al.l (120001 ) and we find their study detects 50% of the total AGB population, and 
their calculated C/M ratio is based on 18% of this populati on. Further, the number of IR 
detected AGB stars we find f or WLM in I Jackson et al.l (120071 ) is in excellent agreement with 



the recent near-IR study of IValcheva et al.l (120071 ). who find a C/M ratio of 0.58. 



We find the total MLRs from AGB stars to be 0.2-1.0x10"^ Mq yr~^, depending on 
the assumed effective temperatures and chemical composition of the mass-loss winds of the 
AGB stars. The distribution of bolometric luminosities and MLRs of individual AGB stars is 



in excellent agr eement with those observed in the Magellanic Clouds (Ivan Loon et al.lll999 
van Loonll2006l ). 
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We show that the dominant contribution (85%) to the thermal near-IR luminosity of 
IC 1613 is from the AGBs and RSGs. 

Finally, the properties of the AGB populations in IC 1613 (i.e., the broadband and 
narrowband optical completeness fractions and MLRs) are nearly identical to those found in 



Jackson et al.l (120071 ) for WLM. This is no doubt due to their similar galaxy characteristics 



and evolutionary histories, and is likely representative of star-forming dis. 

We sincerely thank Serge Demers for providing his narrowband photometry of WLM and 
IC 1613 and for extremely valuable comments regarding narrowband optical completeness 
statistics. We also thank the anonymous referee for their careful reading of the manuscript 
and their valuable comments, particularly regarding the possible detection of extrinsic carbon 
stars. This work is based in part on observations made with the Spitzer Space Telescope, 
which is operated by the Jet Propulsion Laboratory, California Institute of Technology under 
NASA contract 1407. Support for this work was provided by NASA through Contract 
Numbers 1256406 and 1215746 issued by JPL/Caltech to the University of Minnesota. This 
research has made use of NASA's Astrophysics Data System Bibliographic Services and 
the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet Propulsion 
Laboratory, California Institute of Technology, under contract with the National Aeronautics 
and Space Administration. 
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Table 1. Basic Properties of IC 1613 



Quantity 


Value 


Reference 


Right Ascension, a (2000) 


01 04 46.4 


1 


Declination, 5(2000) 


+02 08 45.9 


1 


Heliocentric velocity, V0 (km s^^) 


-232 km s"^ 


1 


Distance, D (Mpc) 


0.73± 0.02 


2 


Morphological Type 


Ir V 


3 


12 + log(0/H) 


7.62 


4 


Total H I mass (M©) 


6.5x10^ 


5 


Inclination angle (degrees) 


38 


1 


Position angle (degrees) 


58 


1 


Rotational velocity (km s~^) 


25 


1 


Conversion factor (pc/arcmin) 


212 


2 



Reference s. — fll iLake fc SkillmanI fll989h. r2llDolDhiii et al. 
J200lh. (3) Ivan den Berghl J2000ah. (4) llee et all J2003h. f5) 
Volders fc HogbomI (1196 ll ). 
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Table 2. Detection Statistics 


Total Point Source Detections in All Wavelengths 


Pilfpr 

J. llLCi 


iXFl 1 TTl n^^T* 


UVJ til V CtilLl 1 


7966 


3.6 


■^91 1 


4.5 /xm 


3162 


5.8 fim 


886 


n //TTi 


618 


Dpfpptinrm in All Fonr TRAC! h^^nrlt; Rnt Not V j^nrl T 






IN LliliUd 


Total 


183 


RSG (Above the AGB hmit) 


3 


AGB/RSG (Above the TRGB, below the AGB limit) 


178 


Below the TRGB 


2 


3.6 //m Point Source Flux (Total Flux = 72.4 mJy) 




Object type 


Fraction 


Brighter than the TRGB 


85% 


Fainter than the TRGB 


15% 


Optical Detection Fractions of IR Identified AGBs 




Filter 


Fraction 


V and I 


57% 


Detected but misidentified in V and I 


11% 


Narrow-band optical 


18% 
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Fig. 1. — Digitized Sky Survey image of IC 1613, 14'xl4' in size. The region of sky with 
coverage in all four IRAC bands is overlaid in black, and the small triangular region with no 
optical Fand / coverage is shown in the southeast. North is up and East is left. 
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Fig. 2.— IRAC 3.6, 4.5, 5.8, and 8.0 //m images of IC 1613. North is up and East is left. 
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Fig. 3. — Absolute 1 versus V—I color magnitude diagram of IC 1613 using data from 



Udalski et al.l (120011 ) . The diagram is segregated into regions containing (a) blue objects, 
(b) AGB stars, (c) red supergiants, and (d) sub-TRGB red giants to aid in determining 
the stellar types of stars in the infrared color- magnitude diagram (Figure H]). The TRGB is 
located at M/ = —4, in the center of the gap separating regions b and d. 
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Fig. 4. — 3.6 yum absolute magnitude versus [3.6] — [4.5] color magnitude diagram of IC 1613. 
The left panel shows all of the po int sources de t ected at 3.6 and 4.5 /im. The right panel 
shows objects also detected in the lUdalski et al.l (l200ll ) V and / data, segregated according 
to Figure [3] with sub-TRGB red giants as red squares, AGB stars as green triangles, RSGs 
as black stars, and blue objects as blue circles. Photometric 1-a errors, averaged over 1 
mag bins, are shown on the right side of the left panel. The AGB limit, TRGB, and a 
vector showing 10 magnitudes of visual extinction are shown in the right panel. See §4] for 
discussion. 
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Fig. 5. — The 3.6 /iin luminosity functions for all of the stars in IC 1613 that were detected 
at 3.6 and 4.5 fim{solid line) and for those detected at 3.6 and 4.5 fim and optically classified 
as sub-TRGB red giants {dashed line). The value we adopt for the 3.6 fxm TRGB is shown 
at M3.6=— 6.2 (see f|l]for discussion). 
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Fig. 6. — The 8.0 /xm absolute magnitude versus [3.6] — [8.0] color magnitude diagram of 
IC 1613. Photometric 1-a errors, averaged over 1 mag bins are shown on the right. 
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Fig. 7. — Positions of stars according to their optical classification, with the exception of the 
bottom right panel, which shows the positions of AGB stars detected in the IR. 
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Fig. 8 . — Fraction of objects more luminous than the TRGB at 3.6 /im that are also detected 
in the lUdalski et al.l (120011 ) Fand / photometry as a function of [3.6] — [4.5] color {left panel) 
and as a function of wind optical depth assuming AGB stars with wind composition of 85% 
AMC + 15% SiC and effective temperatures of 2650 K {top axis) and 3600 K {bottom axis). 
See §5.31 for a detailed discussion of mass-loss modeling. The error bars were created by 
taking the square root of the number of optically detected AGBs and dividing by the total 
number of AGBs seen in at 3.6 and 4.5 /im and therefore only represent the degree to which 
small number statistics might affect the overall trend. 
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Fig. 9. — Absolute 3.6 /im versus [3.6] — [4.5] color-magnitude diagram. The small black 
points are our IRAC d etections and the la rge red circles are our measurements of the car- 
bon stars identified by I Albert et al.l (120001) • Note in particular the four carbon stars with 
[3.6] — [4.5] > 0.2, which are likely losing significant mass. The right panel is magnified to 
help distinguish the points. 
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Fig. 10. — The mass-loss rate versus bolometric magnitude for all objects brighter than 
the the TRGB at 3.6 ^m. The mass-loss rates assume a gas-to-dust ratio of 4.8x10^^, a 
wind composition of 85% AMC + 15% SiC and effective temperatures of Tg// = 2650 K 
(filled circles) and Tg// = 3600 K {open circles), usin g the mass l oss prescription of Ivan Loon 



( l2000l ) (see §5.3p . The bottom dashed line is the lJural (Il984l ) single-scattering mass-loss 
l imit, and the top soli d line is the empirical maximum mass-loss limit observed in the LMC 
( van Loon et al.lll999l ). The AGB limit is shown as a vertical line at Mhoi = —7.1. 



